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Introduction: Ray-tracing (RT) of Lunar Digital 
Elevation Models (DEM)'s is perfomred to virtually 
derive the degree of radiation incident to terrain as a 
function of time, orbital and ephemeris constraints [1- 
4], This process is an integral modeling process in 
lunar polar research and exploration due to the present 
paucity of terrain information at the poles and mission 
planning activities for the anticipated spring 2009 
launch of the Lunar Reconnaissance Orb iter (LRO). 
As part of the Lunar Exploration Neutron Detector 
(LEND) and Lunar Crater Observation and Sensing 
Satellite (LCROSS) preparations RT methods are used 
to estimate the critical conditions presented by the 
combined effects of high latitude, terrain and the 
moons low obliquity [5-7], These factors yield low 
incident solar illumination and subsequently extreme 
thermal, and radiation conditions. 

The presented research uses RT methods both for 
radiation transport modeling in space and regolith re- 
lated research as well as to derive permanently shad- 
owed regions (PSR)'s in high latitude topographic min- 
ima, e.g craters. These regions are of scientific and 
human exploration interest due to the near constant 
low temperatures in PSR's, inferred to be < 100° K. 
Hydrogen is thought to have accumulated in PSR’s 
through the combined effects of periodic cometary 
bombardment and'or solar wind processes, and the 
extreme cold which minimizes hydrogen sublimation 
[8-9], RT methods are also of use in surface position 
optimization for future illumination dependent on- 
surface resources e.g. power and communications 
equipment [10-11]. 

However, ray-tracing methods are computationally 
slow due to the inherent vector transport process and 
the subsequent search phase required to evaluate, iden- 
tify, and interact with surfaces [12]. The issue is exac- 
erbated if the DEM point surface is converted to a 
plate model due to the large number of floating point 
operations (48) required to evaluate each plate facet 
[13-14]. As a result, the computational complexity 
incurred in evaluating direct illumination without op- 
timization constraints increases exponentially as a 
ftmction of the product of the number of surfaces n in 
the model evaluated against all other surfaces during 
search x ~ \n\-l, floating point operations per surface 
evaluation /; number of vectors projected per pixel jv|, 
0{nsjv). In this case, the exponential time-complexity 


increases effectively limit the scale or resolution of the 
model size irt|. In this research, we effectively mini- 
mize the search space size [xl and floating point opera- 
tions I/I, by using vector-point operations, thereby sig- 
nificantly improving the algorithm run-time. 

Methods: Significant computational performance 
gains are obtained by not triangulating the DEM sur- 
face points into plates. Instead, vector-point opera- 
tions facilitate fast evaluation for direct illumination 
and secondary illumination processes, lines (1-9). In 
this approach DEM surface points are warped to a 
sphere (Moon 1734.8 km radius) illustrated in Figure 
1. Given point pj as an arbitrary pixel in image DEM 
for which illumination conditions are required, each 
point must be evaluated for the possibility of shading. 
The Sun is oriented at Cartesian position sun. V| de- 
fines the vector between the sun and pj (1). sun to pi 
distance = tiisti (2), (3) nonnalizes v,. (4) defines a 

linear set dvj of distance scalars from sun as interpo- 
lants along v; at interval, a (km), dv; range = [r/ij/ - 
range, disi^. jdVjj = range ! a. (5) A discrete set of 
Cartesian positions pos e V; are detenuined by project- 
ing from sun along Vj at distances defined by set dry 
(6) npos defines the number of positions in pos to be 
evaluated for shading of p,. (7-11) evaluates each 

point in pos as to its elevation by rounding its (x,y) 
coordinates to common coordinates that index the ele- 
vation scalar in v and DEM. In conditions where the 
elevation z value in the DEM is higher than the z value 
in V, a ‘shade’ condition is identified. 

These steps provide a nearly direct access method- 
ology that significantly minimizes the search space for 
evaluating illumination for each p,, evaluating only- 
possible occluding positions along vector, V|. Starting 
evaluation of pos at dist,, and working towards the sun 
along Vj, facilitates a run time minimization of search- 
ing dVj by monitoring local elevation conditions and 
shaded status to terminate the loop. A consequence of 
this approach is that actual intersection points of v, 
with the surface are not determined. These may be 
approximated using interpolation with sun distances to 
nearest neighbor points. 

Secondary vector projection processes require sur- 
face intersection points and surface normal informa- 
tion. A given vector-surface intersection point poSj is 
obtained by proceeding through each element in pos 
starting at p, to determine the closest shading point p„. 
Surface normals at each point are pre-calcu!aled for 



each model by taking the vector cross product between 
adjacent points, These are averaged to yield DEM 
normals for each point. 

h v, = p, -sun 

2. dist, = [Vj| 

3. Vi = Vj/t/w7i 

4. dv, = interpol(t/jA7„ range, a) 

5. pos = sun r V*dVi 

6 . npos - range/a 

1. f ory - 1 , «/3iw do 

8. r = round(pos(0J)) 

9, y - roTj.nd(pos( I J)) 

10. if DEM(xj.') gt pos(2j) then 

‘Shade 

11, endfor 



Figure 1. Example of single vector direct solar illumi- 
nation. Illumination process is optimized by minimiz- 
ing the number of floating point operations for surface 
evaluations as well as the ray-tracing search space. 
Multiple vector projections are also supported 

Results: Figure 2 illustrates both permanently shad- 
owed and illuminated regions of the Clementine North 
Polar DEM containing 80-90 deg north latitude [15- 
16]. Image resolution (606x606) pixels = 1km*/ pixel 
during summer where the north pole is inclined !.54° 
towards the sun, yielding maximum north polar illumi- 
nation. In this approach, the sun position was rotated 
in 1° increments around the DEM for 360 degrees 
yielding 360 images. Without further optimization, 
each point’s direct illumination is evaluated against all 
points « for shading incurring |«[‘ = 367236 evalua- 
tions to determine the maps direct shading for each 
image. Only a .spaced points along the incident solar 
vector V are evaluated vt'ith maximum required search 
map corner to corner range - 808 points. By integrat- 
ing through the image stack in each pixel, degree of 
illumination is determined. Highliglited PSR's gener- 


ally occur in the bottoms of craters below the southern 
rim, area=1661 km^ (red). Note, linear and blocky 
structures in the DEM are false positives incurred by 
elevation mis-registrations in the overlay of the 
Clementine source imagery during DEM formulation. 
Updates to these products will increase the resolution 
and accuracy of OEM's, modifying the stated PSR and 
illuminated areas. 



Figure 2. North Polar Permanently Shadowed Regions 
derived from Clementine DEM using ray-tracing meth- 
ods for 360, P rotation increments during North Polar 


summer when polar inclination = 1.54°, yielding maxi 


mum illumination. Demonstration used the described 


enhanced ray-tracing methods to derive illumination. 
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